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ABSTRACT 

We consider the transport of angular momentum in accretion disks that are low 
mass, in the sense that the gravitational forces produced by the material in these 
disks has a negligible effect on disk dynamics. There is no established consensus 
on how this transport takes place. We note that for phenomenological reasons the 
traditional a model is probably not a good description of real disks. Here we briefly 
review a few of the more promising models. These include models in which angular 
momentum transport is driven by shocks, and by magnetic field instabilities. The 
latter is more promising, but requires a dynamo. We note that the direction of 
angular momentum transport due to convection in a conducting disk is not known as 
competing mechanisms are at work. We briefly discuss a number of possible dynamo 
mechanisms, and their problems. We then give a detailed exposition of the internal 
wave driven dynamo model, in which internal waves excited at large radii drive an 
a — SI dynamo. The azimuthal magnetic field produced in this way is unstable to a 
magnetic shearing instability (MSI), which drives nearly isotropic turbulent eddys 
with typical fluid velocities of ~ Va , where Va is the local Alfven speed. The scale 
of the eddys is ~ Va/SI, where fi is the local rotational frequency. This turbulence 
leads to a saturation of the dynamo when Va ~ (H/r) 2 ' 3 c s , where H is the half- 
thickness of the disk, c s is the local sound speed, and r is the radial coordinate. This 
gives rise to an effective dimensionless viscosity coefficient ~ (iJ/r) 4 / 3 and vertical 
and radial diffusion coefficients which are ~ (H/r) 4 / 3 Hc s . The resulting vertical 
diffusion of entropy will have a substantial effect on detailed models of vertical 
structure in accretion disks. Viscous and thermal instabilities of very hot disks, 
those dominated by radiation pressure and electron scattering, are substantially 
moderated in this model. We note that the MSI largely suppresses the Parker 
instability in accretion disks. 



1. Introduction 

The study of astrophysical disks has, at its center, a curious and enduring 
puzzle, i.e. the radial transport of angular momentum. Throughout this book such 
transport is assumed to take place at a rate much higher than one would obtain 
by appealing to well understood microscopic processes. The evident success of such 
models, in spite of their ad hoc basis, suggests that some kind of collective process 



(or processes) is responsible. Here we will describe recent work aimed at producing 
a model of these processes. What follows will not be completely objective. There 
exists no well-defined consensus in this field and no model is free from serious 
objections. This paper describes what we think are the most promising approaches 
to this problem. Other workers will have different opinions. A popular, although 
somewhat dated, reference for this topic is Pringle 1981, who quite forthrightly 
declared that very little was known about disk viscosity at that time. We believe, 
perhaps optimistically, that times have changed. Nevertheless, much of the following 
discussion will be qualitative, which is indicative of the large uncertainties that 
remain. 

We begin by defining the problem we wish to solve in the most narrow pos- 
sible terms. First, we are concerned with thin accretion disks, i.e. those for which 
the disk thickness H is much less than the radius r. This is not because we con- 
sider thick disks to be rare or uninteresting, but because (H/r) is a natural ordering 
parameter. When it is large then the problem becomes significantly more compli- 
cated. Second, we assume that the local rotational frequency, Q(r), is essentially 
Keplerian, i.e. Q 2 = GM*/r 3 . If the disk is thin, then this is equivalent to assuming 
that M dlsk < M*. Third, we will assume that our disks have negligible self-gravity. In 
other words, their vertical gravity is approximately -z£!(r) 2 , where z is the distance 
to the midplane. This implies that M dlsk < (H/r)M*. A disk that violates this con- 
dition will suffer from local gravitational instabilities which are outside the scope 
of this paper. We note that as a consequence of these three conditions c s ~ HQ, 
where c s is the sound speed at the disk midplane. Fourth, we ignore any externally 
imposed magnetic field. Fifth, we assume that the disk is perfectly conducting. 
These conditions need not be met everywhere within a disk. Rather we expect that 
we can confine our attention to those regions within a disk where these conditions 
are met. Viewed in this light these are not particularly onerous conditions and will 
apply to the majority of visible systems. The chief exceptions are protostellar disks, 
which are likely to be poorly conducting. We will ignore their existence. 

The traditional treatment of accretion disks relies on the introduction of a 
phenomenological viscosity given by v = ac s H (Shakura and Sunyaev 1973). This ex- 
pression was originally motivated by the expectation that the high Reynolds number 
of accretion disk flows would lead to strong, approximately isotropic, turbulence. 
Since H and c s are the only natural scales for length and velocity this implies an 
a of order unity, the exact value of which is only weakly dependent on the specific 
nature of the disk in question. Since then it has become clear that there are no 
purely hydrodynamic linear instabilities in accretion disk flow. Nevertheless there 
still may be some justification for this approach. Processes which are purely local 
in the sense that they do not rely on input from distant regions of the disk will 
tend to reduce to this formalism since H and c s remain the only natural scales. 
This includes local MHD processes, since the only 'natural' scale for the magnetic 
pressure is the pressure of the conducting fluid. Quite apart from the question of 
whether or not a models will turn out to be justified due to some local instability, 
they have been extremely useful in constructing phenomenological models of disks. 



This is due, in part, to the fact that the description of a stationary, optically thick 
disk is independent of the precise nature of the disk 'viscosity'. Still, real disks 
are not necessarily either stationary or optically thick and it should be possible 
to constrain models of angular momentum transport by appealing to observational 
work. In general, nonlocal models will not be equivalent to a models, but it is still 
convenient to quote results in terms of an equivalent "a" for a stationary disk. 

It is worth noting that it is not clear that a purely local theory that is really 
equivalent to an a model can be made to work at all. A number of attempts have 
been made to model dwarf novae outbursts as the result of a thermal instability in 
accretion disks around white dwarfs and neutron stars (Mineshige and Osaki, 1983; 
Meyer and Meyer-Hofmeister, 1984; Meyer, 1984; Smak, 1984; Mineshige and Osaki, 
1985; Lin, Papaloizou, and Faulkner, 1985; Cannizzo, Wheeler, and Polidan, 1986). 
Similar work has also been done on X-ray transients (Huang and Wheeler, 1989; 
Mineshige and Wheeler, 1989). Such models can produce quite successful results, 
but only if the value of a can be made to rise significantly as the disk temperature 
goes up. Moreover Wood and Mineshige (1989) have argued that during quiescence 
the emission profiles of the disks of Z Cha and OY Car imply that a decreases at 
small radii as r °- 3 ~ 0A . These results can be summarized as suggesting that a scales 
as H/r to some power of order unity. It is not clear how a model based on local 
processes can produce this kind of scaling. 

A useful way to look at angular momentum transport processes is to start 
by examining the available linear MHD modes, and see what roles they might play. 
We do this in section 2 of this paper, at the same time reviewing several candidate 
processes. In section 3 we will describe the basic theory of magnetic dynamos in 
disks and discuss some previous attempts to construct a self-consistent disk dynamo. 
In section 4 we return to the focus of our current work, the internal wave driven 
dynamo and show how to derive the critical scaling arguments that describe its 
operation. In the final section we will discuss the prospects for progress in this 
contentious field. 

2. Local Modes 

Our discussion of the local modes within a disk will be based on a number of 
simplified situations, which will suffice to illustrate the critical physical processes. 
We begin by making a methodological point. There are a variety of equally valid 
ways to examine modes in a system which is undergoing some zeroth order motion. 
One common method used in this context is to describe disturbances in terms of 
some radial plane wave, i.e. f(r,6) oc exp[i(k r r + m6)]. The vertical structure may be 
included explicitly, giving a nontrivial dependence on z, or, for some applications, 
its effects may be ignored and the perturbation may be assumed to be proportional 
to exp(ik z z). In either case the underlying differential rotation of the disk will give 
the perturbation a complicated dependence on time, including a linear increase in 
k r with time. We will prefer to describe perturbations as rigidly rotating patterns, 
with a simple dependence on time, but a complicated radial structure. In general 



we feel that this makes the evolution of waves and instabilities more transparent. 

We start by considering local adiabatic hydrodynamic disturbances in a disk. 
In this case 'local' means that the typical radial length scales associated with the 
perturbations are much less than the radius of the disk. We will therefore neglect all 
radial derivatives of the background fluid properties and the geometry. There is one 
important radial derivative that survives (and explains why the radial eigenfunctions 
are not pure plane waves). This is because d t always appears in the combination 
d t + fld e = dt + imQ. This can be understood physically as being due to the fact 
that the physically important frequency is not the one measured by an observer 
in an inertial frame, but rather the frequency of the perturbation measured by an 
observer moving with the unperturbed fluid. This is ui = w + mil. A mode with a 
given frequency will have properties, determined by u>, that change over a typical 
distance of ru>/(mfi). 

The general linear equation governing such disturbances is 
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The use of \ has no particular significance except that it simplifies the relevant 
equations. N is the Brunt- Vaisala frequency. It vanishes at the midplane of the disk 
and will be large (compared to fi) and positive in the disk atmosphere. In between 
it may have any value and will be imaginary in convectively unstable regions. The 
associated velocity perturbations are 
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and the associated density perturbations are 
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A qualitative sense of the solutions for these equations may be obtained by replacing 
the first two terms in eq. (la) with -k 2 x- 



Equations (1) and (2) describe linear perturbations, but in order to estimate 
the consequences of various linear processes we need to have some way of estimating 
the nonlinear saturation of these modes. In general this can be quite difficult, and 
requires a detailed examination of the mode-mode couplings. In the particular case 
of rapid, small scale strong turbulence we will adopt the usual heuristic approach 
of assuming that the effect on other modes can be modeled by a turbulent diffusion 
tensor of the form = (viVjT corr ). If the turbulence is isotropic then we can ap- 
proximate this by DSij where D is the mean square velocity of the turbulence times 
the eddy turnover time. There is an additional dissipative term which is unique 
to shearing environments (Vishniac and Diamond 1989). Radial diffusion coupled 
to differential rotation will smear wave motions over an azimuthal wavelength at a 
rate of ^ g 

TsD~((j) 2 n 2 D rr ^ . (3) 

We will refer to this process as the shear decorrelation rate. It can be important in 
setting limits on the importance of high to modes in accretion disks. 

The angular momentum flux, C z is given by pr{v r vg). On physical grounds we 
expect that for axisymmetric waves this is entirely due to mixing. Since turbulent 
mixing of a conserved quantity tends to minimize its gradient, and since specific 
angular momentum increases outward in a thin accretion disk, the net effect will 
be to move angular momentum inward. In fact, from eqs. (2a) and (2b) we see that 
for ™ = 

(^) = 2| _;f^ 2|2 i^xi 2 ^), (4) 

i.e. growing modes have C z < as expected. 

The energy flux associated with disturbances in the disk can be shown (cf. 
Vishniac and Diamond 1989) to have the form 

F = (SPv r ) + prfl{v r vg) + £{pv r ), (5) 

where £ is the energy density. The last term describes the advection of energy by 
mass transport and is of no particular concern. The first one is the transport of 
energy expression that one would obtain in the absence of a strong zeroth order 
flow. It gives the transport of energy measured by a local observer, i.e. terms 
that go as the square of the perturbation variables. This term points in the same 
direction as the wave group velocity, as one would expect since it represents the 
transport of a positive quantity by the waves. The second term is just C z Vi. Waves 
carrying angular momentum will slow down, or speed up, the fluid while they pass, 
and this gives rise to an energy contribution proportional to the product of rQ. and 
v g . Since C z can have either sign the total energy flux may, in general, have either 
sign as well, i.e. a wave propagating inward may have a positive total energy flux. 
To linear order both F and L z are preserved as a wave propagates. This has the 
consequence that (5Pv r ) is not. It is particularly important to note that if {SPv r ) 
and £ z have opposite signs then as a wave propagates inward fl will increase and so 
will (SPv r ). Under most circumstances {v 2 ) will as well. 



It is also interesting to note that the condition that C z > (or in other words, 
that the waves carry angular momentum outward, thereby driving the evolution of 
the accretion disk) is that u? increase in the direction of wave propagation. Inward 
propagating waves are then particularly helpful since they have (5Pv r ) < 0, C z > 0, 
and, since d r fl < and u> = w + mfl, will tend to grow in amplitude as they move 
inward. This should promote nonlinear dissipation and the deposition of negative 
angular momentum into the disk fluid. This is the basis for theories of angular 
momentum transport using compressional modes (Shu 1976) and internal waves 
(Vishniac and Diamond 1989). We will now examine each in detail. 

2.1. Sound Waves 

When perturbations have large frequencies, in the sense that Co 2 > ft 2 , then 
neglecting vertical structure we can write 

w 2 » c 2 s k 2 + WKB terms, (6) 

i.e. the usual one for sound waves. The WKB terms represent the secular effect of 
the shearing environment on the wave amplitude, a point we have already mentioned 
in connection with the separate conservation of energy flux and angular momentum 
flux. A detailed examination of these waves shows there is a lower limit to lj 2 for 
a given k z ^ and m. As a sound wave reaches this limit it will undergo reflection. 
No corresponding upper limit appears. It follows that a sound wave has, at most, 
one surface for radial reflection. It appears from this, and from the preceding 
discussion, that one might construct a model of angular momentum transport based 
on these waves in the following manner. If there is some source for sound waves 
at the outer disk boundary, that produces waves that can travel through the disk 
without reflection and that have C z > then nonlinear interactions between the disk 
and sound waves may allow the deposition of negative angular momentum and the 
inward transport of mass. This is, in essence, the proposal of Shu (1976) who pointed 
out that the accretion stream creates a shock wave when it hits the disk, and the 
resulting spiral shock may penetrate deeply into the disk. Of course, the properties 
of such a shock depend to a large degree on nonlinear terms neglected in eq. (1) 
nevertheless it can still be used as a crude basis for understanding the results of 
more detailed calculations. Following Shu's original suggestion, analytical work by 
Donner (1979), Spruit (1987) and Larson (1990) has shown that the dimensionless 
value of a obtained in this way scales as (H/r) 3 ^ 2 . This is significantly larger than 
the estimate one would obtain from examining the angular momentum flux carried 
by a sound wave and reflects the ability of coherent, nonlinear effects to enhance the 
efficiency of angular momentum transport. A number of numerical studies (Sawada, 
Matsuda, and Hachisu, 1986; Sawada, Matsuda, Inoue and Hachisu, 1987; Spruit, 
Matsuda, Inoue and Sawada 1987; Rozyczka and Spruit 1989; Matsuda, Sekino, 
Shima, Sawada, and Spruit 1990) have confirmed that spiral shocks result from 
both tidal effects and the impact stream, and that the shock waves penetrate to 
very small radii. 



In spite of the many attractive features of this model, several problems re- 
main. First, almost all this work has been done in two dimensions. The only 
exception has been Sawada and Matsuda (1992) who ran a 3 dimensional simu- 
lation of the formation of a disk from an accretion stream far enough to see the 
formation of the spiral shocks, but not nearly far enough to show a steady state 
solution. In a two dimensional disk N 2 = and c s is a constant. We see from eq. 
(1) that the resulting linear solutions are almost trivial and neglect several physi- 
cally important features. Moreover they have the unpleasant feature that since \ is 
not a function of z the velocity perturbations « c s x-P (1/7)_1 diverge in an isothermal 
atmosphere. This is the result of refraction, a point examined by Lin, Papaloizou 
and Savonije (1990). They showed that generic compressional disturbances inside 
an accretion disk will tend to quickly reach the disk atmosphere and undergo rapid 
nonlinear dissipation. Realistic solutions to eq. (1) will also show this effect, but 
have a compensating feature. As mentioned above, all waves undergo reflection as 
N 2 increases past u>. This means that a global sound wave with an azimuthal wave 
number m will undergo reflection at ~ m scale heights. For m small this will keep 
the wave from dissipating in the disk atmosphere. Unfortunately, this advantage 
will be lost as the wave steepens into a shock and its characteristic time scale falls 
to zero. We can only conclude that vertical confinement of the wave energy seems 
to be a serious problem. 

Second, this process is relatively ineffective. The scaling of a with (H/r) 3/2 is 
consistent with most observational constraints, but the coefficient in front is quite 
small, — 0.01 (Spruit 1987). The result is a level of angular momentum transport 
which cannot be easily reconciled with our current understanding of dwarf novae. 
Some numerical work suggests that this low value is not correct (Matsuda, Sekino, 
Shima, Sawada, and Spruit 1990), but more recent work by Rozyczka and Spruit 
(1992) is consistent with the earlier analytic estimates. 

Third, Livio and Spruit (1991) have examined the outburst intervals of dwarf 
novae as a function of mass ratio as a test of the scaling of a with (H/r) and concluded 
that a gentler dependence might be preferred. Evidently there are few grounds for 
believing that spiral shocks mediate angular momentum transport in dwarf novae 
systems. If we conclude that this process does not dominate in dwarf novae then 
it becomes hard to understand why it should dominate in other, less well observed, 
systems. 

2.2. Internal Waves 

The second set of hydrodynamic modes are internal waves, analogous to deep 
ocean waves. For small m they have a local dispersion relation 



This is identical to the dispersion relation for deep ocean waves, except that we 
have replaced the coriolis parameter with Q 2 . In practice, the difference is that in 
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the ocean we normally consider the coriolis parameter to be small, whereas £1 2 is a 
large parameter for internal waves in an accretion disk. We note that in the limit 
N 2 = k z = this dispersion relation collapses. Purely two dimensional treatments 
of disk structure will not include these modes. For three dimensional disks this 
dispersion relation implies that internal waves are confined between a radius of 
reflection, where uj 2 — > fi 2 and k 2 — > and a radius of resonant absorption where 
lj 2 — > and k 2 — > oo. Everywhere in between internal waves are confined vertically to 
regions where u? > N 2 . As a wave approaches reflection the group velocity implied 
by eq. (6) is proportional to the square root of the distance to the reflection point. 
Consequently reflection takes a finite amount of time, as expected. Near reflection, 
a wave occupies a broad slice of the disk midplane, although it still cannot propagate 
into regions where N 2 > tt 2 . At the other extreme, as a wave approaches resonance 
its group velocity is proportional to the square of the distance to the resonant radius. 
Consequently the wave energy piles up and undergoes nonlinear dissipation. This 
is exacerbated by the fact that as a wave approaches resonance it is confined to 
an increasingly thin zone around the midplane. We can conclude from this that in 
general these waves are transient events that are confined to annuli with a width of 
order r/m. Regardless of starting conditions they eventually enter a state of rapid 
nonlinear dissipation, close to their corotation radius. Stationary modes based on 
eqs. (1) or (6) are unphysical in that they contain an infinite amount of energy at 
their corotation radius. 

This generic picture does not apply to modes with m = 0,±l. The former 
case is of limited interest, since such modes contain no angular momentum (and are 
also incapable of driving a dynamo). In the latter case, it is possible to arrange the 
parameters of the waves so that their resonant points lie outside the disk and their 
reflection point is at r = 0. To be more specific, if we choose waves with u less than 
Q, at the outer edge of the disk and require m and u to have opposite signs, then 
q = u> - Q fulfills these conditions. Such waves will be global modes of the disk. This 
allows us to return to the picture we discussed earlier, in which angular momentum is 
transported outward by inward propagating waves whose linear amplification leads 
to strong nonlinear interactions and the deposition of negative angular momentum 
at small radii (cf. Vishniac and Diamond, 1989). It turns out that in conducting 
disks this effect is overwhelmed by the role of these waves in driving a dynamo. Still 
it is worth discussing some of the details of the propagation of these waves, both 
to compare the results to the spiral shock wave model, and to set the stage for a 
discussion of the wave-driven dynamo model. 

Assuming, for the moment, such waves are excited at large amplitude at the 
outer edge of the disk, what happens to them as they move inward? First, they 
will undergo linear amplification at a rate of V group /r. Modes with small m and 
uj/Q of order unity have V group ~ c s (fc z i?) _1 . The amplification will also be driven 
by geometric focusing as H and r decrease together, although this will increase 
the result only by a factory of order unity. In a steady state this increase will 
be balanced by nonlinear dissipation. The nonlinear dissipation of such waves is 
complicated. The radial group velocity has a strong dependence on q so we expect 



that incoherent nonlinear interactions will be important. This suggests a generic 
nonlinear dissipation rate of 
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which is consistent with work on deep ocean waves (McComas and Bretherton 
1977; Miiller et al. 1986) There may be important coherent effects as well, e.g. the 
formation of solitons, but little is known about this possibility or its consequences. 
Equating the linear amplification rate with r^ corr then for u> some fraction of il we 
find that 

(v 2 ) ~ *c 2 s {k z H)-\ (9) 

We see that the wave spectrum will be dominated by modes with k z H ~ 1. The 
corresponding angular momentum flux is 

C z ~ ^c 3 s . (10) 

Assuming that nonlinear dissipation deposits the angular momentum in the waves 
into the disk at a rate of C z T^ corr /V group gives an 'alpha' equivalent to 

a~(H/r) 2 (11) 

There is a major assumption contained in these estimates, that the nonlinear 
interactions between the various modes lead to a steady state with (u> 2 /Q 2 ) a constant 
less than one everywhere inside the disk. If this fails then u> — > as the waves 
go to small radii and the group velocity and amplitude will both decline sharply. 
Simulations of the evolution of the internal wave spectrum based on the weakly 
resonant approximation and assuming random phases supports this assumption 
(Huang 1992), but more realistic calculations are still desirable. 

Is it fair to assume that internal waves of this strength are excited at the 
outer edge of the disk? It is true that waves of smaller strength will grow to 
saturation given enough time, i.e. a large enough number of radial e-foldings, but 
this growth goes as r to a modest power and will be insufficient if the waves are very 
weak at the outer edge of the disk. At this point the excitation of internal waves 
is not understood, and there are a number of unresolved issues associated with this 
problem. A non-self-gravitating disk in a binary system has two possible sources 
of negative angular momentum, tidal interactions and the accretion stream. The 
former is generally supposed to be the ultimate sink of excess angular momentum 
from the disk, and will involve a loss of angular momentum at a rate of ~ Mr 2 n. The 
latter is due to the addition of material with a specific angular momentum far below 
the average value at the outer edge of the disk, so that the stream will constitute a 
local sink of angular momentum with a total rate comparable to the tidal loss rate. 
An ingoing stream of \m\ = 1 internal waves with an amplitude of (v 2 ) ~ (H/r)c 2 s (see 
below) will carry an angular momentum flux of approximately T,H 2 c 2 . This implies 
that either the accretion stream or tidal interactions need to excite ingoing |m| = 1 



internal waves with an efficiency of ~ (r/H) 2 a. Goodman (1992) has suggested that 
a standing wave pattern of m = -1 waves should result from the action of the m = 2 
component of the tidal forces, but this may be sensitive to the inner and outer 
boundary conditions of the disk. Lubow (1992) has pointed out that tidal forces 
will also excite vertical motion within the disk, which corresponds to the excitation 
of the minimal k z mode of internal waves, but no quantitative results are available. 
In two dimensional models the impact stream excites modes with u = uj binary and 
k z = (i.e. x not a function of z), but in real life a substantial component of 
the forcing will have \ a function of z and drive internal waves. This question 
has not been properly addressed since previous treatments of the interaction of 
the disk with the accretion stream and/or the tidal forces from the secondary star 
have been almost entirely two dimensional. Even granting the existence of one, or 
more, effective excitation mechanisms there is also the issue of how quickly high 
amplitude, coherent waves excited near the outer edge of the disk evolve into the 
expected low amplitude, incoherent waves we expect at smaller radii. In sum, there 
seems to be plausible grounds to suspect that these waves are excited, but almost 
nothing is actually known about the process. We note that the issue here goes well 
beyond the existence of the waves. If tidal forces are ineffective at exciting them, 
but the accretion stream provides a strong excitation mechanism, then temporary 
interruptions in the accretion stream ought to result in a rapid cooling of the disk. 
Such interruptions seem to occur frequently in cataclysmic variables just above the 
period gap (Robinson, Barker, Cochran, Cochran and Nather 1981), and do result 
in a rapid descent into an extremely quiescent state. If we consider an isolated 
disk then the situation is even murkier. Adams, Ruden, and Shu (1989) and Shu, 
Tremaine, Adams, and Ruden (1990) have shown that massive disks (i.e. disks with 
a total mass comparable to the central object) have a fairly robust m = 1 instability 
centered on a narrow annulus at the outer edge of the disk. More recently Noh, 
Vishniac, and Cochran (1992) showed that a variant of this instability can persist 
down to very low mass disks, for which M dlsk ~ (H/r)M*, i.e. disks for which our 
(defining) constraints on mass are marginally violated. Unfortunately the coupling 
between this instability and ingoing internal waves has not been explored. In any 
case, isolated disks for which M dlsk < (i?/r)M* even at the outer edge would seem to 
be bereft of any exciting mechanism. Fortunately, the notion that such disks are 
completely inactive seems consistent with all the available observational evidence. 

How does this estimate compare to the one derived from the theory of spiral 
shock waves? The scaling law presented here is steeper, which may be even more 
of a problem for dwarf novae models (cf. Livio and Spruit). On the other hand, 
the coefficient, although uncertain, ought to be of order unity, suggesting that this 
effect will dominate for disks with (H/r) greater than about 1CT 4 . The vertical 
structure of the disk poses no particular problem for internal waves, in contrast to 
the difficulties and uncertainties involved with sound waves and shocks. On the 
other hand, we have already noted that the excitation of internal waves is much 
less well understood. Finally, the complicated nature of the nonlinear evolution of 
the internal waves simply makes it more difficult to provide reasonable analytical 



estimates of the coefficients of the scaling laws in this case. 



2.3. Convection 

It is immediately apparent from eq. (7) that when N 2 < an instability is 
lurking in the dispersion relation for internal waves, i.e. standard convection in a 
shearing environment. The presence of shear poses no particular obstacle to this 
instability (although it does limit it in some interesting ways) and so we would 
expect convective turbulence to result. This has been invoked many times as the 
basis for a theory of angular momentum transport (Cameron, 1969; Paczynski, 
1976; Vila, 1978; Lin and Papaloizou, 1980; Lin, 1981; Smak, 1982; Cabot et al, 
1987a,b). 

Although it is clear that convective mixing ought to transport angular mo- 
mentum, there are some difficulties with using it as the basis for a general model 
of angular momentum transport. For example, convection will continue only if the 
entropy of the gas continues to decrease away from the midplane. The entropy 
gradient is, in turn, a function of opacity and the vertical profile of heat deposition. 
In this model, heating is due to the dissipation of orbital energy, which is, in turn, 
due to convective mixing. There is no obvious reason why this should form a self- 
consistent cycle. Also, optically thin, nearly isothermal disks seem to exist, and in 
such disks a convective cycle will not start, and if started will quickly die away. 

Perhaps the central question concerning convection, however, is that of the 
direction of angular momentum transport. For example, Goodman and Ryu (1992) 
have asserted that convection invariably transports angular momentum inward in 
accretion disks, on the basis of a shearing sheet calculation in which k r — ► oo as 
t oo. This is clearly a crucial issue, as only outward transport will facilitate the 
accretion process. We have already noted that for axisymmetric modes angular 
momentum is transported inward in an accretion disk. The ambiguity in the sign of 
C z is a natural consequence of the fact that in an accretion disk, the driving entropy 
gradient is vertical, not radial. Hence, only vertical entropy mixing follows from the 
standard lore concerning the relation between gradients and the instabilities they 
cause. The intrinsic difficulty in determining the angular momentum flux, as well 
as its relation to the convective cell dynamics, is nicely illustrated by considering 
Rayleigh-Bernard convection in a thin disk in solid body rotation, with |7V | 2 of order 
Q, 2 . In this case the convective growth rate is 

2 _ {k 2 . + (m/r) 2 )\N\ 2 ~k 2 z m 2 

k 2 + k 2 + (m/r) 2 ' 1 ' 

where we have replaced uj 2 with -r 2 since we are dealing with an instability whose 
frequency is no longer a function of r. The quantity 40 2 = 2£lr~ 1 <9 r (r 2 0) is the 
Rayleigh discriminant (equal to O 2 for a Keplerian disk). Equation (12) expresses 
the well-known relation between convective cell structure and instability enforced 
by the Taylor-Proudman theorem for a rapidly rotating fluid. Instability requires 
that the vertical cell scale associated with vertical vortex line bending not be less 



than the planar cell scale, i.e. k\ < k 2 or (m/r) 2 . Note that cells with k 2 > k 2 z > (m/r) 2 
are shaped like toroidal rolls with an ellipsoidal cross section in the r-z plane, while 
cells with (m/r) 2 > k 2 > k 2 resemble anisotropic Proudman Pillars, with an elongated 
ellipsoidal cross section lying in the r-9 plane. The impact of this eddy anisotropy 
on angular momentum transport is in turn illustrated by the quasi-linear expression 
for the angular momentum flux, which is: 



The first two terms are primarily due to coriolis force effects, while the third is 
related to pressure gradient induced transport. The resulting angular momentum 
transport is clearly very sensitive to eddy anisotropy, in that for k r < (m/r) (corre- 
sponding v r > vg for incompressible flow with k z small) transport is outward, while 
for (m/r) < k r (corresponding to v r < v e ) transport is inward. Isotropic convection 
cells make little or no contribution to L z . We note that the third term is, in general, 
quite significant, although unless the symmetry between r and 6 is broken, contri- 
butions from modes with opposite helicity will cancel. We see that the sign of C z is 
determined by the detailed structure of the eddy wave-number spectrum, which is 
difficult to predict. 

In a thin disk differential rotation distorts finite k z modes into locally unstable 
spiral modes. We expect that the linear eigenf unctions of eq. (1) for this case will 
be similar to the non-rotating case sufficiently close to the corotation point that 
< r, i.e. with a scale length of ~ Yr/(mVi). At larger distances the instability 
will look like outgoing internal waves and the amplitude will fall off rapidly. The 
net angular momentum flux is 



We note that the symmetry between (m/r) and k r is explicitly broken in this ex- 
pression. Moreover, since r and 9 are no longer symmetric, there is no reason to 
expect that the third term will tend to sum to zero. Given that the peak of \x\ 2 is 
confined within an annulus of size Tr/m£l we can see that the second term is always 
negligible. In this sense, convection in accretion disks does seem to favor inward 
transport of angular momentum. Unfortunately, when m/r is large, so that the 
radial wavelength of the instability fills the convective annulus, then the third term 
will will be comparable to the first. Since k. r does not become arbitrarily large as 
t — > oo in the normal mode analysis, the exact sign of C z for a given mode can only be 
determined by explicit calculation. Moreover, since large m modes will be favored 
by phase space arguments, and not more subject to dissipation unless m > rk z , the 
general direction of angular momentum transport due to convection in accretion 
disks must be regarded as uncertain. It is interesting to note that this question is 
related to the question of whether or not differential rotation inhibits convective 
instabilities. On energetic grounds we expect that if C z < 0, so that differential 




(13) 



C z = pr( 



oj 2 -n 2 \ 2 



(nr (2 (^) 2 |x| 2 - i|3„ x | 2 ) + (n 2 M a )i^x*a*)>. (u) 



rotation acts as a sink of energy in convective regions, then convective instabilities 
will be inhibited in Keplerian disks. 

Finally, convection will not be terribly effective as a source for a unless it is 
very strong. We note from eq. (7) that coriolis forces limit the radial wavelength 
of a convective cell to less than ~ fcj 1 r/0, where r = iN. For weak convection the 
convective cells will be tall and thin. Approximate incompressibility implies that 
v r ~ (k z /k r )v z and on dimensional grounds we expect v z ~ v$ ~ fe~ 1 r. These estimates 
are self-consistent in the sense that using them we obtain D z - fc~ 2 r and D r - k~ 2 T 
so that the turbulent damping rate is comparable to r. These estimates also assume 
that (m/r) is at least not very much larger than k z . Similar arguments suggest that 
very large m modes will tend to be relatively ineffective in driving transport. On 
the other hand, it is well to remember that this discussion is implicitly quasi-linear, 
and real convective turbulence shows features that are not easily understood in 
such simple terms. As a best estimate we conclude that a <~ (v r v g )/c 2 ~ (r/fi) 3 in 
agreement with Ruden, Papaloizou and Lin (1988) (although their discussion was 
based on axisymmetric modes for which a is necessarily negative). We stress again 
that the actual sign of a in this case is uncertain. 

We note that Korycansky (1992) has suggested that an infinitesimal amount 
of dissipation will eventually suppress any convective mode, albeit after an extended 
period of growth. This is an artifact of relying on a plane wave approximation 
for the radial structure of the convective modes, which includes the unfortunate 
consequence that k r oc t and the dissipative damping rate k 2 D — > oo as t — > oo. 

In view of these problems we tend toward the view that while convection 
may have some effect on disk structure, it is probably limited to a subsidiary role. 
A general model for a must come from some other mechanism. 

2-4- The Magnetic Shearing Instability 

Inasmuch as ionized disks are extremely well conducting systems it seems 
reasonable to explore the possibility that they contain some magnetic field. Such a 
magnetic field has been suggested as the source of angular momentum transport in 
disks by a number of authors (e.g. Lynden-Bell 1969; Shakura and Sunyaev 1973; 
Eardley and Lightman 1975; Ichimaru 1977; Coroniti 1981). Due to the presence of 
differential rotation any radial field component will lead directly to a rapid increase 
in B g . It follows that we should expect the magnetic field to be predominantly 
azimuthal. Perturbations of a highly conducting disk with an azimuthal field will 
have a fairly complicated dispersion relation. Rather than present the generalization 
of eq. (1), which is not particularly enlightening in any case, we will take a simplified 
case. If we assume incompressibility and N 2 = the dispersion relation becomes (cf. 
Vishniac and Diamond, 1992) 
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The other hydrodynamic perturbation variables are 
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and 



where V A = (B 2 /4ttp) 1 / 2 is the Alfven speed. The magnetic field perturbations have 
the form 

b = B e —v + tB e — r v r 9. (16d) 

A somewhat more straightforward dispersion relation results if we assume that we 
are in the limit of large radial wavenumbers, i.e. k r > mfl/(rui) so that we can 
approximate the radial eigenfunctions as plane waves. Then the RHS of eq. 12a 
can be neglected and the surviving terms can be written in the form 

+ - uj) + Co 2 {-Su 2 A uo 2 + l^llltu 2 A n 2 ) - Aujw\ - 0, (17a) 

where 

J\ = k ln 2 . (176) 

In the limit where u? A < uj and m/r < k z eq. (17) has a root with u? w 3(2(m/rk z ) 2 - 
i.e. an instability. Some discussion of this case using a radial plane wave 
WKB analysis can be found in Balbus and Hawley (1992). 

This instability was discovered by Velikhov (1959) and described at length 
by Chandrasekhar (1961) in the context of a vertical magnetic field threading an 
incompressible Couette flow. We will refer to it here as the Magnetic Shearing 
Instability (MSI). Its novel feature is that it depends only on d r 9, < 0, rather than 
an unstable gradient in the specific angular momentum, and necessarily acts to 
transport angular momentum in the direction of decreasing angular speed. It is also 
present in the limit of an arbitrarily weak magnetic field, although its transport 
effects become negligible as B — >• 0. These novel features follow from the way that the 
distorted magnetic field lines serve to tie together fluid elements at different radii. 
The more slowly moving fluid at larger radii is accelerated by fluid elements lying 
on the same field line but at smaller radii, thereby transferring angular momentum 



outward. In spite of the relatively early date of this work its importance for the 
dynamics of accretion disks was not realized until recently (Balbus and Hawley 
1991; Hawley and Balbus 1991). Subsequent work on the nonlinear evolution of 
this instability (Hawley and Balbus 1992a; Zhang, Diamond, and Vishniac 1992) 
has been largely confined to a consideration of the unstable axisymmetric modes of 
a vertical field. This turns out to be of limited use in considering the role of the 
MSI in a realistic accretion disk. 

What are the critical wavelengths for MSI? We note from eq. (16) that the 
instability exists for uu\ < uj, in which case the growth rate is of order cu A , but 
is suppressed for u>\ > u>j. This implies that it has a maximum growth rate of 
order Q, which is achieved only when k z is not much less than |fc| and (m/r)VA ~ 
Apparently the fastest growing modes have azimuthal wavelengths of Va/SI and 
radial and vertical wavelengths which are comparable to each other, and no larger 
than the azimuthal wavelength. More slowly growing modes can have much larger 
azimuthal and vertical wavelengths. We note that they cannot have larger radial 
wavelengths. The validity of the plane wave approximation rests on the condition 
that the radial wavelength be less than the characteristic distance for changes in u>. 
For MSI acting on an azimuthal magnetic field this implies 

d r ui (m/r)n Q , . 

LU UJ A V A 

In other words, MSI modes for an azimuthal field have a (narrow) characteristic 
radial scale. This seems at first to be counterintuitive. After all, one could always 
define initial conditions with a large radial scale. Equation (18) follows from the fact 
that after one characteristic growth time differential rotation will have automatically 
reduced that initial scale, no matter how large, to ~ V A /£l. This argument does not 
apply to the MSI modes of a vertical field, which have a growth rate proportional to 
k z VA instead. In that case, the vertical wavelength is limited only by the thickness 
of the disk, and the radial wavelength is limited to a similar value by secondary 
instabilities (Zhang, Diamond, and Vishniac 1992). 

What is the angular momentum transport associated with the MSI? We have 

£z = pr(v r v e ) - r— . (19) 

Using eqs. (15) and (16) we can show that 

_^ oj a \ 1 iCu 2 \d r v r \ 2 



We note that an unstable mode will be localized around the corotation point, where 
Cu is purely imaginary, and will have an amplitude that drops sharply to zero as 



the real part of Q becomes important. In this region all the terms comprising C z 
are positive for a growing mode, except for the first, whose sign is ambiguous. 
Fortunately, the lower limit on k r imposed by differential rotation implies that 
this term is smaller than the others. The magnetic shearing instability transports 
angular momentum outward. 

Although a detailed evaluation of the nonlinear effects of the MSI has not 
been completed for the physically interesting case of an azimuthal magnetic field, 
some qualitative comments are in order. We can begin by considering the modes 
for which (m/r) ~ k z ~ k r ~ SI/Va and which consequently grow at rate ~ fl. The 
instability will tend to induce velocity perturbations of order V A . This speed is, of 
course, what one would expect from a process driven by magnetic forces, but we also 
note that this is sufficient to alter the local value of the shearing thereby eliminating 
the driving force behind the instability. Since the wavenumbers are all comparable 
the local turbulence induced on these small scales will be approximately isotropic. 
The turbulent diffusion coefficients will all be of order V%/Sl. Similarly, we can see 
from eq. (19) that the angular momentum flux will be — V\ for an effective 'a' of 
{Va/cs) 2 . There remains the question of whether or not the fastest growing modes 
are actually the ones that dominate the transport terms. More slowly growing 
modes with longer characteristic scales would seem to have an advantage. In the 
case of an azimuthal field this does not appear to be the case. Since all modes have 
comparable radial wavelengths, nonlinear interactions between modes appear in the 
force equations with an associated time scale of k 2 r D r ~ (fl/V A ) 2 V^/il ~ O. For all but 
the fastest growing modes this term will totally dominate the dynamics of motion, 
destroying the long azimuthal wavelength modes through turbulent diffusion. We 
are led to a picture in which a large scale magnetic field coexists with a turbulent, 
disordered field of comparable amplitude on scales of V A /tt whose features change 
at a rate ~ f2. 

This argument does not apply to the MSI of a vertical field, in which the 
slowest modes have k r ~ k z ~ l/H and a turbulent diffusion rate of H~ 2 V%/Vt ~ 
(Va/H)(V a /c s ). In this case, the turbulent diffusion rate is smaller than the lin- 
ear growth rate by a factor of Va/c s and the slow linear modes are unaffected by 
small scale turbulence. Eventually the azimuthal field will build up to the point 
where its MSI can suppress the vertical field modes. This implies a saturation field 
strength such that V\ ~ Va z c s , with typical turbulent eddy scales of ~ (VAzH/fl) 1 ^ 2 . 
By contrast, in two dimensional numerical models the azimuthal field has no ef- 
fect on the instability and grows until it reaches V A ~ c s (Hawley and Balbus 1991; 
Zhang, Diamond and Vishniac 1992). In either case, the resultant angular mo- 
mentum transport is equivalent to an a of ~ Va z /c s (Zhang, Diamond and Vishniac 
1992). Preliminary results from three dimensional simulations are consistent with 
our prediction of a low saturation value for V A (Hawley and Balbus 1992b), but 
have insufficient dynamic range to provide a conclusive test. (We note in passing 
that these authors take an opposing view of the interpretation of their numerical 
experiments.) 

The MSI seems exceptionally well suited to angular momentum transport in 



an accretion disk. It is a robust instability which moves angular momentum in the 
right direction. It even allows us to estimate the dimensionless viscosity. The work 
presented above, however, is not a complete predictive theory. Unless we have a 
way of understanding the equilibrium strength of the magnetic field in the disk we 
have merely substituted one free parameter, a, for another, B\/P. To go further, 
we need to understand the accretion disk dynamo. 

2.5. Magnetic Buoyancy 

One important effect associated with vertical motions in a magnetized fluid 
is magnetic buoyancy, or the Parker instability. Although quite general, the details 
of this effect depend on the detailed vertical profile of the fluid and its derivation 
is rather involved. Here we will simply quote some relevant results and refer the 
reader to Parker (1979) for a detailed treatment of this process. 

The basic physical driving mechanism for this instability is the buoyant force 
associated with the fact that gas with an associated magnetic field can maintain a 
higher pressure, for a given density, then unmagnetized gas. In a magnetized fluid 
rising field lines can shed mass if other sections of the same field lines are sinking, so 
that gas frozen to the field lines can slide off the rising magnetic bubbles. Depending 
on the detailed vertical profile of the magnetic field and the gas this instability will 
lead to a bending mode of the magnetic field lines. This instability can actually 
be enhanced in the presence of high thermal conductivity. The vertical wavelength 
of the unstable modes is on the order of a scale height. The growth rate can be 
as high as Va/H for modes with a wavelength parallel to the field lines of order 
H. The instability is insensitive to the transverse wavelength. It follows that this 
instability can persist in accretion disks if we restrict ourselves to sufficiently short 
radial wavelengths, i.e. k r > Q/(i2T) ~ fl/V A . This was confirmed in a detailed 
calculation by Shu (1974). 

The usual picture of this instability is that it leads to flux expulsion at a 
rate of Va/H as rising magnetic loops reconnect in the disk atmosphere. Vishniac 
and Diamond (1992) pointed out that this may be extremely misleading. Although 
magnetic buoyancy persists in the linear theory of MHD perturbations in accretion 
disks, the usual dimensional estimates of vertical mixing and magnetic flux ejec- 
tion depend on the saturation state of this process. However, in accretion disks a 
significant complication arises. The minimum radial wavenumber associated with 
magnetic buoyancy is, as noted above, ~ CI/Va- This wavenumber is virtually the 
same as the one associated with the MSI so we can expect a close interaction be- 
tween the two sets of modes. Worse, the Parker instability is relatively slow for 
« c s . If we assume that the action of the MSI on the Parker instability can 
be represented by a kind of turbulent diffusion, then the associated diffusive term 
gives a damping rate of ~ kf(V^/fl) — ft. Physically, what is happening is that the 
MSI modes are causing the rising segments of field lines to trade momentum with 
nearby sinking segments on a time scale much shorter than the growth time for the 
instability. Of course, the buoyancy of the field cannot be entirely suppressed, but 



the associated growth rate is reduced to r ~ (T 2 Parker /(k 2 D)) ~ (V A /c s ) 2 tt. Assuming 
efficient reconnection, this rate is the rate at which the Parker instability can eject 
magnetic flux from the disk. This revised estimate is not only much smaller than 
one obtains by ignoring MSI, it is also comparable to the rate at which turbulent 
diffusion (due to MSI) will transport flux. Evidently the magnetic field buoyancy 
does not introduce any new time scales. 

We also note that shear decorrelation effects will also suppress the Parker 
instability. From eq. (3) we get for (m/r) > H^ 1 



This is still faster than the magnetic buoyancy rate, although not as fast as the 
radial mixing rate. 

This does not amount to a complete suppression of the Parker instability in 
accretion disks. It will still compete with turbulent diffusion in determining the 
vertical gradient of the magnetic field. Also, we would still expect the emergence 
of magnetized bubbles to contribute to heating the disk atmosphere. The strong 
interaction between the MSI and the Parker instability does imply that magnetic 
buoyancy does not act to limit the growth of magnetic fields in disks. Moreover, 
the loss of spatial and temporal coherence of the modes of the Parker instability 
makes it a poor candidate for the driving mechanism in a magnetic dynamo, a role 
that has been previously suggested by Stella and Rosner (1984), and more recently 
by Tout and Pringle (1992). We will return to this point in the next section. 

2.6. Boundary Instabilities 

Although no local hydrodynamic instabilities are evident in eq. (1) it turns 
out that one can induce instabilities by attaching the accretion disk to a rigid 
radial boundary. This has been shown by numerous authors (Drury 1980, 1985; Pa- 
paloizou and Pringle 1984, 1985; Goldreich and Narayan 1985; Narayan, Goldreich 
and Goodman 1987; Glatzel 1988) in various ways. The instability arises from the 
ability of reflecting radial boundary conditions to cause over-reflection of the im- 
pinging waves (Goldreich and Narayan 1985), and vanishes for absorbing boundaries 
(Blaes 1987). If a natural reflection boundary for the waves exists inside the disk, 
then this process can be repeated until they shock. It is possible that this plays 
a role in the dynamics of boundary layers in accretion disks, but it is extremely 
unlikely that this is important in the bulk of the disk. Numerical studies ( Glatzel 
1986; Goldreich, Goodman and Narayan 1986; Hanawa 1988a,b) have shown the 
effects of this instability are concentrated near the boundary, while calculations by 
Kaisig (1989) have demonstrated that this process saturates in a series of weak 
shocks that have only a small effect on angular momentum transport. 
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2.7. Finite Amplitude Instabilities 



In spite of the stability of accretion disks to infinitesimal hydrodynamic dis- 
turbances (see also Dubrulle and Knobloch 1992) there remains the possibility that 
some finite amplitude disturbance might manage to perpetuate itself, or generate 
equivalent motions nearby. This possibility has been explored through work which 
examines the response of the disk to a local defect in the flow (Lerner and Knobloch 
1988; Dubrulle and Zahn 1991) and also through studies of simple turbulence clo- 
sure schemes (Iroshnikov 1980; Dubrulle and Valdettaro 1992; Dubrulle 1992). In 
the first case one can see that the defect induces instabilities in the surrounding 
flow. The second approach has yielded the suggestion that an inverse turbulent 
cascade may arise in accretion disks which prevents the destruction of large scale 
defects in the flow and maintains a moderately strong turbulent flow. It is difficult 
to evaluate the future prospects for this work. In spite of the work done to date the 
persistence of defects in the flow cannot be rigorously justified. On the other hand, 
it is difficult to rule out such effects. Current three dimensional hydrodynamic work 
has not confirmed this model, but problems with dynamic range and artificial vis- 
cosity still limit the ability of computers to confront these models directly. If this 
work is correct, then it falls into the category of a local model of angular momentum 
transport for which a has a definite and universal value. 

3. Disk Dynamos 

The generation of large scale magnetic fields in highly conducting fluids, 
like stars or accretion disks, is not particularly well understood, although it seems 
clear from stellar magnetic reversals that this process must be possible. In what 
follows we will adopt the formalism of mean-field dynamo theory, while mentioning 
unresolved issues associated with this approach. Reviews of mean-field dynamo 
theory can be found in Moffatt (1978), Parker (1979), and Zel'dovich, Ruzmaikin 
and Sokoloff (1983). The basic idea in this approach is that the magnetic field 
can be divided into a large scale, slowly evolving component, B, and a small scale, 
rapidly changing component, b. The local velocity field interacting with B, through 
the induction equation, will induce some b. This same velocity field acting on b will 
produce a change in B. A large scale, fast dynamo results if there exists a mode 
with B increasing. In order to obtain the mean field dynamo equations one iterates 
through the induction equation and averages over small distances and times. We 
obtain 



where e^k is antisymmetric under the interchange of indices, repeated indices implies 
summation and 
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The Dij tensor is the usual diffusion tensor. The a„ tensor is known as the helicity 
tensor and its trace is the helicity of the velocity field. The derivation of eq. (21) 
assumes an incompressible velocity field, which is a reasonable approximation in 
the context of the motions proposed to drive an accretion disk dynamo. 

Equation (21) is based on the assumption that the back reaction from the 
small scale field is negligible and that the large scale field evolution time scale is 
long compared to the relevant time scales for b and v. In some circumstances b can 
have a dramatic effect on v, and consequently on the growth of B, and this can lead 
to a saturation of the dynamo process. This can usually be taken into account by 
modifying the definition of a tj . A more serious problem is that the second term 
in eq. (21a) has a diffusive role, analogous to the (neglected) effects of ohmic 
dissipation (Parker 1971; Vainshtein and Ruzmaikin 1971, 1972). This actually 
represents a diffusion of the likely location of large scale magnetic field lines, but in 
order to use this 'turbulent diffusion' as a substitute for ohmic dissipation we need 
to believe that small scale reconnection of field lines actually removes local sharp 
gradients in the field rather than merely moving them around in an unpredictable 
manner. Recently Vainstain and Cattaneo (1992) have pointed out that if this 
reconnection is assumed to take place on scales characteristic of ohmic dissipation 
then the amount of small scale magnetic field energy required by the reconnection 
of large scale field lines is prohibitively large. The back reaction from such fields 
would lead to a saturation of the dynamo process when the large scale field strength 
was almost completely negligible. The Kraichnan spectrum (Kraichnan 1965) for 
MHD turbulence predicts little energy on such scales and suggests that large scale 
field lines will be relatively unbent by a strongly turbulent cascade if the magnetic 
field energy and fluid kinetic energy are roughly equal on large scales. Since large 
scale fast dynamos do appear to exist this poses something of a paradox, whose 
resolution is not entirely clear. One possibility for the Sun is that relatively large 
scale loops can be created and transported by convective flows, so that reconnection 
events consist of the mutual annihilation of large scale, oppositely oriented field lines 
(Vainshtein, Parker, and Rosner 1993). This should also work in accretion disks 
where the MSI will constantly create loops of size V A /Ct and differential rotation 
should help increase field gradients. The reader should bear in mind that this 
plausible suggestion is not supported, either in the Sun or in accretion disks, by 
detailed calculations. 

In a differentially rotating disk eq. (21) must be modified to take into account 
the background flow of material. This has the effect of adding the term rd r QB r 6 to 
the right hand side. In an accretion disk this term will dominate the growth of Bg. 
This leads to a cycle in which fluid motions, acting through the a tj term generate 
B r from B g . The differential shearing of the disk then closes the loop, and the 
dissipative terms provide the damping. This is referred to as an a - Q dynamo. 
The total dynamo growth rate will be less than O (since the generation of B r from 
Bg will tend to be inefficient). This in turn implies that the dominant large scale 
field component will be Bg. For a thin accretion disk the only important spatial 
gradient in B (or in and D l} ) will be in the z direction. This allows us to produce 



a simplified dynamo equation of the form 



d t B r = -d z (a g0 B ) + d z (D zz d z B r ), 
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and 



d t B e 



-flB r + d z (D zz d z B e ), 



(226) 



with age and D zz as defined in eqs. (21). Note that we have dropped the third term 
in eq. (21a). As long as B is largely azimuthal this term will be very small. A more 
exact treatment would include the a er B r term in eq. (22a), but such a term cannot, 
by itself, drive a dynamo since a 2 dimensional dynamo is impossible (Zel'dovich 
1956). The requirement that VB = implies that B z is of order (H/r)B r for the 
internally generated field in a disk, which makes it dynamically negligible. 

From eq. (22) we can see that if a growing mode exists its growth rate will 
be approximately 



Of course, this estimate is only realistic if we can neglect the damping terms. The 
existence of a growing mode requires that 



where a gg and D zz are suitably averaged over the vertical structure of the disk. 
Equation (24) assumes that the vertical scale height of all the relevant physical 
quantities is H. More realistically, one has to solve the appropriate vertical structure 
equations with a given model of the small scale velocities. 

So far we have not addressed the question of the velocity field that drives 
age, although we have mentioned several candidates in previous sections. We will 
go through the list of candidate processes starting with those least likely to drive 
a dynamo successfully. We will defer a discussion of the wave driven dynamo until 
the next section. First on this list is the Parker instability. The unstable modes 
of magnetic buoyancy will have k z ~ (m/r) ~ H -1 , and k r ~ CI/Va- Turbulent 

mixing due to MSI implies a coherence time for these modes of The vertical 
acceleration will be of order V%/H so the vertical velocity after one coherence time 
will be V|/c s . Since the flow is approximately incompressible this implies v r ~ 
{k z /k r )v z or v r - (V A /c s ) 3 c s . It follows that 



The minimum reasonable value for D zz is the one given by the MSI, or V|/fi. Com- 
paring eqs. (24) and (25) we can immediately see that magnetic buoyancy cannot 
drive a successful dynamo when H <c r. 

The next candidate mechanism is the MSI. There is almost no work on the 
ability of the azimuthal MSI to drive a dynamo, although Zhang et al. (1992) 
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have shown that the vertical MSI can do so. This result is interesting from the 
point of view of supplying a seed field for some other dynamo process, but by itself 
appears to generate an azimuthal field whose final strength depends on the strength 
of the applied external field. As far as the azimuthal MSI goes we can only make 
some tentative and qualitative comments. We start by noting that eqs. (15) and 
(16) contain neither vertical structure nor local buoyancy. Although this should be 
adequate for modes with rapid growth rates and small vertical wavelengths, it leaves 
us with a set of relationships that give no indication of any loss of local symmetry 
that would make a e e nonzero. If we explicitly calculate aee we obtain 



where T corr is the correlation time for the mode, which in the saturated state should 
be roughly the inverse of the growth rate. This expression should be zero (after 
averaging over a set of similar modes) for two reasons. First, the real contributions 
come from terms that are odd functions of the distance to the corotation point, 
and will integrate to zero after radial averaging. Second, there is no obvious reason 
why one sign of k z should be preferred. This argument has the loophole that it 
is implicitly weakly nonlinear, whereas the saturated state of the instability will 
leave these modes in a strongly nonlinear state. Nevertheless, the basic lack of a 
means of generating the necessary asymmetries remains. We might instead appeal 
to modes with k z H ~ 1, which have the possibility of significant contributions from 
buoyancy effects and vertical structure terms, but we can see from eq. (17) that the 
maximum growth rate on these scales is — (k z /k)Q <~ Va/H. Such modes will have 
k z - (ra/r) ~ H^ 1 and k r - 0,/Va- They will be suppressed for the same reason that 
magnetic buoyancy modes are suppressed, and by an identical argument can be 
shown to be incapable of driving a dynamo. Of course, if either of these processes 
could be shown to work in spite of these difficulties then one would have a theory 
of angular momentum transport which would be a true a model, in the sense that a 
would be a constant of order unity describing a purely local process of dissipation. 

It turns out that convection is a particularly interesting candidate for driv- 
ing a dynamo, a point first raised by Galeev, Rosner and Vaiana (1979). When 
convection is weak, with a growth rate r c « J] we have seen that it is characterized 
by k z ~ (m/r) ~ H^ 1 and k r ~ (il/r c )_ff _1 with v z and v e of order HT C and v r smaller 
by a factor of T c /Q. When the MSI mixes a convective cell faster than r c then 
convection will be disrupted and, like the Parker instability, convection becomes a 
poor prospect for the source of a successful dynamo. On the other hand, when 



then convection will proceed unimpeded. In this case, vertical mixing occurs at 
a rate of r c so that D zz - H 2 T C . The correlation time of the convective flow is 
approximately r^ 1 so 




(26) 




(27) 




(28) 



We see from these results and from eq. (23) that convection is marginally capable 
of driving a dynamo. Even if it is unable to drive a large scale dynamo it is clear 
that the convective cells will tend to wind up the magnetic field lines of a disordered 
field. The local field strength will consequently grow until the MSI associated with 
it is strong enough to disrupt the convective cells. When eq. (27) is replaced by 
an equality then the angular momentum transfer associated with convection will be 
approximately equal to that associated with MSI. The difference is that the MSI 
will clearly induce an outward transfer of angular momentum, unlike purely hydro- 
dynamic convection whose effect is uncertain. It seems possible that sufficiently 
strong convection in a conducting medium transports angular momentum outward 
with a dimensionless 'viscosity' a ~ (T c /il) 3 . In addition, turbulent diffusion and 
magnetic buoyancy would tend to spread the magnetic field generated within a 
zone of strong convection until such fields were found for all z at any radius where 
convection occurs. Such fields would, of course, be subject to the MSI and promote 
outward angular momentum transport. 

4. The Internal Wave Driven Dynamo 

At this point we have a good candidate for the local driving mechanism 
behind angular momentum transport (the MSI), but the best dynamo candidate 
discussed above is convection, which will not always exist, or necessarily sustain 
itself self-consistently, and cannot be guaranteed to drive angular momentum out- 
ward instead of inward. Clearly we need an alternative. One attractive possibility 
is that the internal wave field discussed in section 2.2. might drive a dynamo. A 
similar point was originally raised by Pudritz (1981a,b) who envisioned the waves 
as arising spontaneously from some local instability in the disk. This process does 
not seem to occur, but the global excitation of such modes can stand in its place. 
Our results here are taken from Vishniac, Jin and Diamond (1990) and differ from 
Pudritz's treatment not only in the source and amplitude of the wave field, but also 
in the details of the derivation of a ee from such waves. 

4-1. Internal Waves as a Source of Helicity 

We have already argued that the dominant internal wave modes have k z ~ H^ 1 
and u largely real. Starting from eq. (2) we find that 

k —id 2 

aee * s2{ w^)W~m dzX * x)Tcorr - (29) 

We have assumed m small, but not not necessarily ±1 since nonlinear interactions 
will generate other values of m. This expression will be nonzero for any given mode, 
but unless the sign of mk r is fixed the contributions from different modes will cancel 
out. Fortunately the sign of mk r is fixed by the expectation that the wave field 
will be dominated by modes that are undergoing linear amplification. We can see 
from eq. (7) that internal waves have a radial group velocity with the opposite 



sign of their phase velocity, i.e. sign[V gr ] = sign[w/k r ]. The requirement that internal 
waves are undergoing amplification is that they be headed for reflection, i.e. that 
w 2 increase in the direction of propagation, implying that sign[V gr du) 2 } > 0. This in 
turn tells us that 

3m 

Luk r dCu 2 = ttuj 2 k r > 0. (30) 

r 

We conclude that mk r < for modes undergoing linear amplification. The existence 
of strong differential rotation in the disk fixes the sign of mk r of the dominant wave 
modes, which in turns determines the sign of their contribution to agg. 

It is useful to rewrite eq. (29) in terms of the mean square velocity contained 
in modes with a given m and u. Unless u> is very close to fl we have (v 2 ) ~ (Sxk z /u) 2 
so that 

agg ~ (v 2 ) ( — J k r —HT corr , (30) 

where we have neglected iV 2 . Including it would confine the waves, and the helicity, 
to the region around the midplane where uj 2 > N 2 , but within that region properly 
accounting for iV 2 will introduce a correction of order unity. Vishniac and Diamond 
(1992) have shown that including the confinement of the waves does not affect the 
vertically averaged angular momentum transport significantly so we will ignore this 
point in what follows. For realistic vertical structure agg will be an odd function of 
the distance to the midplane of the disk. For strongly turbulent flows T corr is the 
eddy turnover time. For waves it is 3?(iw) _1 . This is zero for waves with exactly real 
frequencies. Purely periodic motions do not result in any net transport. For finite 
amplitude waves, however, nonlinear interactions between the waves give rise to a 
small imaginary frequency, of order r^ orr . Consequently, r corr « (<i' 2 T (;ecorr )~ 1 . We 
conclude that 

aee - (v 2 K e l rr (I') kr^H. (31) 

If we imagine that the energy from the fundamental, global wave mode is passed 
through a weakly turbulent cascade down to small scales, then within that cascade 
we expect (u 2 )r^ orr integrated over a logarithmic interval to be roughly constant. 
To the extent that further linear amplification can take place within this cascade 
( v2 ) T de 1 corr m ight even increase. Consequently eq. (31) suggests that the helicity will 
be dominated by modes with large m and k r . Of course, this argument fails if the 
sign of mk r is not preserved in the cascade, but preliminary calculations (Huang 
1992) suggest that it is preserved. 

4-2. The Internal Wave Spectrum 

Given that a gg is likely to be dominated by small scale waves, generated from 
the nonlinear dissipation of the global modes, we need some understanding of the 
course of nonlinear dissipation. Here we appeal to a simple model for the cascade, 
discussed in Vishniac, Jin, and Diamond (1990) and Vishniac and Diamond (1992). 
Previous work on the dissipation of deep ocean waves (McComas and Bretherton 



1977) has shown that there are three basic processes mediating the nonlinear trans- 
fer of energy in a weakly turbulent cascade: elastic scattering, induced diffusion, 
and parametric subharmonic instability. The first two processes are unimportant in 
this context. The third involves the generation of waves from a wave with roughly 
twice the radial wavenumber and half the frequency of the generated waves. The 
vertical wavenumber changes only slightly in this process (remember that this is a 
continuum model!). The analog for an accretion disk is that the global modes tend 
to generate local modes with k z ~ H^ 1 but steadily decreasing w and increasing k r . 
This process will eventually terminate in the eruption of strong turbulence at some 
suitably small scale. This model makes no firm prediction as to the sign of mk r in 
the cascade or the magnitude of m as a function of k r . Numerical work by Huang 
(1992) based on the weakly resonant approximation (Phillips 1960, 1961) supports 
this basic picture and also suggests that m increases slowly in the cascade while 
the sign of mk r is preserved. In what follows we will take (v 2 )T^ corr as a constant 
in the cascade. This can be justified in detail by comparing T decorr and the time 
required for linear amplification. This implies that (v 2 ) oc ui 1/2 and T decorr oc w 1/2 . 
Since (v 2 ) - (H/r)c 2 for the global modes which lie at the top of the cascade we have 

<A> ~ ?^(£V /2 , (32a) 



H /^ 1/2 



and 

irr,u> ~ 7 (J J ( 326 ) 

We will also assume m 2 oc but our results change only slightly if the increase in 
m is more gradual. We obtain from eqs. (31) and (32) 

h\ 3 /^ 3/2 



UJ Cs (33) 

All of the previous work cited is based on a random phase approximation, so our 
use of it here implies a neglect of soliton formation in the ingoing internal waves. 

In an unmagnetized disk the cascade will extend down to scales where 
(v 2 ) 1 / 2 > Huj. On smaller scales the fluid will show strong, anisotropic turbulence. In 
a magnetized disk these results will be valid only if d> > (m/r)VA and if the turbulence 
due to MSI dissipates the waves more slowly than r^ corr . Below this scale the fluid 
motions will be rapidly dissipated by the small scale turbulence induced by the 
MIS. The former condition is necessary for our treatment of internal waves as purely 
hydrodynamic motions to be valid. It is also easy to satisfy. The latter condition is 
important in providing a natural saturation mechanism for the dynamo. From eq. 
(32b) we can see that it is equivalent to 

^i<" 2 (f)(§) 1/2 > (34«) 



or 



, r x 2/3 /V A \ 4/3 
*>(-) (£) 0. (346) 



As Va increases this limit will come to define the lower end of the weakly turbulent 
wave cascade. Since a gg is dominated by waves at this end of the spectrum the 
increase in Va will have the effect of turning down the dynamo growth rate and 
increasing the turbulent dissipation rate until they balance. 



4-3. Scaling Laws for the Internal Wave Driven Dynamo 

We are now in a position to derive the basic scaling laws that describe the 
internal wave driven dynamo model. When the dynamo is saturated, the dynamo 
growth rate, given in eq. (23), must be balanced by turbulent dissipation. Since 
the MSI induces roughly isotropic turbulence this implies that 

^dynamo ~ ^' $5) 

From eqs. (23) and (33) we see that T dynamo is dominated by the contribution from 
wave modes with the smallest values of u> that still satisfy the inequality given in eq. 
(34b). Combining these results we find that in a stationary disk, with a saturated 
dynamo, 
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and 
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2/9 

n, (36d) 



where ui c is the frequency for which eq. (34b) is an equality. We see that this fre- 
quency is rather insensitive to the shape of the accretion disk, and will be small, 
relative to ft, only for extremely flattened systems. We have also neglected correc- 
tions due to the vertical structure of the disk and its ability to confine waves with 
small lj to regions very close to the midplane. An approximate correction to a due 
to this effect can be found in Vishniac and Diamond (1992). It changes eq. (36c) 
by a factor of only (iZ/Azc) 1 / 6 , where Az c is the height of the region where waves of 
frequency u> c can propagate. For a nearly adiabatic disk this is a factor of ~ 1, but 
even for an isothermal disk this factor is only (H/r)~ 2 / 15 . 

The model of the turbulent spectrum that we have used to arrive at these 
results is somewhat crude. Still, we can see that the exact cutoff frequency and 
dynamo growth rate are not terribly sensitive functions of (H/r) and will not be very 
different if we had assumed m constant (or increasing as w _1 ) within the turbulent 
cascade. Taking a somewhat broader view, we can see that no model of the turbulent 
cascade, even one incorporating soliton effects, will give a outside the range of 



(H/r) 1 ^ 15 . At the upper limit for the exponent we have the contribution due to the 
global modes alone. At the lower limit the MSI turbulence will start to interfere 
with the propagation of those modes, and thereby saturate the entire process. 

At this point it is useful to summarize the basic physics of the internal wave 
driven dynamo. The flowchart below illustrates the way the various elements of 
this model interact to drive mass inward in an accretion disk. Properly speaking 
the scaling laws we have quoted here apply only to a stationary system. For disk 
evolution rates in excess of Y dynamo the arrows on this flowchart must be replaced 
by time dependent equations. For evolution rates in excess of the wave propagation 
rate, (H/r)Cl, deviations from these stationary estimates can be large. 

EXCITATION OF INTERNAL WAVES 

(+ Linear Amplification) 

| (+ Nonlinear Wave Interactions) 
SATURATED WAVE SPECTRUM 
(+ Shearing) 
a - n DYNAMO 

GROWTH OF B e 

(+ Magnetic Shearing Instability) 
SMALL SCALE TURBULENCE 



SATURATION TURBULENT TRUNCATION SUPPRESSION 
OF B e TRANSPORT OF INTERNAL OF PARKER 

WAVE SPECTRUM INSTABILITY 

4-4- Phenomenology of the Internal Wave Driven Dynamo 

We have seen that the internal wave driven dynamo gives a result for the 
angular momentum transport in a stationary disk which is apparently consistent 
with phenomenological models. This is actually not certain since this mechanism is 
non-local and may show significant deviations from a local model with a ~ (H/r) 4/3 . 
At this point not much is known about where such deviations might show up. Here 
we comment on a few phenomenological aspects of this model that can be derived 
from qualitative considerations. 



It is well known that a model accretion disks dominated by electron scattering 
whose main pressure support is due to radiation are subject to thermal and viscous 
instabilities (Pringle, Rees, and Pacholzcyk 1973; Lightman and Eardley 1974). 
Lightman and Eardley also noted that these instabilities can be cured by excluding 
the radiation pressure from the definition of viscosity, but this suggestion has no 
obvious physical basis. It is certainly not a consequence of the internal wave driven 
dynamo model. Fortunately, it turns out that these instabilities are largely, though 
not entirely, suppressed in this model. In order to see this we need to review the 
basis of these instabilities. The rate at which local thermal equilibrium is established 
in a disk, rf 1 , is the local flux divided by the vertical integral of the thermal energy. 
The rate at which a stationary surface density distribution is established, r^ 1 is 
the mass transfer rate divided by Sr 2 . For a stationary disk dominated by electron 
scattering and radiation pressure we have 

M fH s 2 
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— V, (376) 



and 



T 1 ~ ^ ~ afl, (37c) 

where F 1 is the photon flux from the disk and P 7 is the photon pressure near the 
midplane. The first condition is the condition for hydrostatic equilibrium and is 
established at a rate of ~ f2. Consider first a perturbation to the structure of 
the disk in which some annulus gains an extra bit of thermal energy. The heat 
generation rate is just Mfi 2 or a£iJ 2 f2 3 . On thermal time scales £ is a constant, so 
the heating rate grows as aH 2 , where H is proportional to P 7 . On the other hand, 
the cooling rate is 

Fj „ „ (38 ) 



where n is the average mass per electron. This is proportional to H so the heating 
rate grows more quickly than the cooling rate and the disk undergoes a phase of 
runaway heating. Conversely, if we had considered an initial drop in the local value 
of P 7 the disk would, locally, have undergone runaway cooling. This instability can 
be cured if a is proportional to H n , with n < -1. At first glance it would appear 
that the internal wave driven dynamo model, and any model designed to reproduce 
the phenomenology of cataclysmic variable disks, can only exacerbate this problem. 
Moreover, since the annulus of the perturbed region has not explicitly entered our 
discussion, we can expect this difficulty to crop up on all scales larger than H. 

These disks are also viscously unstable, as we can show by considering longer 
time scale variations and assuming that thermal equilibrium is somehow maintained. 



If F 1 ~ MVL 2 then 



H~^. (39) 

Cfl 



Then remembering that M ~ a'EH 2 ^, we see that M oc (Sa) _1 . For a constant a we 
note that overdense regions in the disk have a smaller M and will become increasingly 
dense and cold, while underdense regions in the disk have a larger M and will 
become increasingly diffuse and hot. Again, this instability can be evaded if a oc H n 
with n < — 1. It is not clear what the consequence of these instabilities might 
be. Apparently, they make it impossible to construct a locally sensible theory of 
accretion disks in this limit. 

It turns out that in the internal wave driven dynamo model these instabilities 
are completely suppressed for annuli with widths much smaller than the disk radius. 
We start by considering very short time scales, shorter than r^ J/ 1 namo . From flux 
conservation we see that Bg oc if -1 . Consequently a oc (Va/c s ) 2 oc (i/ 3 S) _1 or a oc H~ 3 . 
We see immediately that this suppresses instabilities on such short time scales. 
Since r t _1 ~ T dynamo , i.e. the critical time scale is just long enough to force us to 
consider the evolution of the dynamo equations, this estimate does not answer the 
general question of the stability of narrow annuli in radiation pressure and electron 
scattering dominated disks. Suppose we consider longer time scales. As long as 
the annulus is narrow compared to r the internal wave energy within the annulus is 
fixed by the condition that the flux supplied from larger r is constant. In an annulus 
that is thicker than its upstream neighbor the waves will spread out and speed up. 
We can express the result in terms of R, the ratio of the mean square perturbed 
velocities contained in |m| = 1 waves relative to the normal saturated value. Since 
Y,(v 2 )c s is constant we have 



Of course, a change in R affects the nonlinear energy transfer rate as R 2 . Tracking 
the effects of a change in R through our previous results we find that w c is unaffected 
but that a is proportional to R. In other words 



Combining eqs. (40) and (41) we see that on thermal time scales a oc if~ 5 / 3 £ _1 . 
Short radial wavelength thermal instabilities are suppressed. On long radial scales 
the internal wave flux can grow (and shrink) to track the saturation value of (u 2 ), 
implying the possibility of thermal instabilities on such scales. Even if these occur 
they do not affect our ability to construct locally sensible models of an accretion disk. 
Such long radial scale instabilities would imply that the disk luminosity would vary 
significantly on time scales — (iJ/r)~ 4 / 3 f2 _1 . These variations could not be arbitrarily 
large since once SH/H was of order unity across the unstable annulus the waves 
would no longer be able to grow, or decay, fast enough to remain near the normal 
saturation limit. 

If we now consider viscous time scales, so that S is no longer held constant, 
we find that M is not a function of £ at all. The internal wave driven dynamo model 




(40) 




(41) 



for radiation pressure and electron scattering dominated disks is neutrally stable 
against short wavelength viscous instabilities. Low levels of absorption and gas 
pressure will stabilize these modes, but on longer radial scales we expect that viscous 
instabilities will be much harder to suppress than thermal instabilities. We therefore 
expect these disks to show luminosity variations on time scales of (H/r)~ w l 3 Q.~ l . 

One point that will be a generic feature of any model that relies on the 
MSI of an azimuthal magnetic field is that the dominant turbulent eddys are nearly 
isotropic. Consequently, the vertical diffusion coefficient, V^/fl, is comparable to the 
radial diffusion coefficient and also to aHc s . This implies that there will always be 
a minimal level of vertical mixing taking place, even in convectively stable accretion 
disks. To see how important this is, we can compare the radiative flux per unit 
area in a stationary disk, Mfl 2 , to the typical vertical heat transport due to mixing. 
The latter is ~ PD Z /L S , where L s is the entropy scale height. The ratio between 
these two transport terms is ~ H/L s , which will usually be of order unity. In other 
words, vertical radiative transport in accretion disks must always compete with 
hydrodynamic mixing effects of the same order. Of course, in a stably stratified 
disk the latter will tend to move entropy inward, implying a greater temperature 
difference between the midplane and photosphere than for purely radiative models. 

Another interesting point about this model is that a fair amount of energy 
is available for coronal heating. The magnetic field energy per unit area in the disk 
is ~ (H/r) 4 / 3 PH. This will diffuse into the disk atmosphere at a rate of ~ (H/r) 4/3 Q 
so we see that the fraction of the energy dissipated locally which is available for 
magnetic heating of the corona is ~ (H/r) 4 / 3 . On the other hand, the internal wave 
energy per unit area is — (H/r)PH. Nonlinear dissipation will dump this energy into 
neighboring modes in phase space at a rate of ~ (H/r)Q. Most internal waves in 
the disk will only couple efficiently to other internal waves, but the global modes 
are not far removed from the lowest frequency sound waves and we can expect 
some fraction of their energy to be transferred to sound waves. Eventually a large 
fraction of the energy in sound waves will end up as shocks in the disk atmosphere 
giving a heating rate of ~ (H/r) 2 PH£l per unit area, or an efficiency for coronal 
heating that will scale as (H/r) 2 / 3 . Although the scaling law for shock heating is 
more favorable, the coefficient is likely to be smaller since it relies on a series of 
inefficient intermediate steps. Unless H/r is very small both processes are likely to 
be important in coronal heating. 

Finally, we note that we have, by hypothesis, excluded external magnetic 
fields from our discussion, but there are certainly circumstances where such fields 
are likely have a large effect on the structure of the disk. In extreme cases, like 
intermediate polars, the magnetic field of the central star may be large enough to 
disrupt the disk entirely inside some critical radius (Ghosh and Lamb 1979). This 
should certainly be the case when the energy density of the imposed field exceeds 
the kinetic energy density of the disk, i.e. B 2 /(Att) > p(rSl) 2 /2. At smaller field 
strengths the presence of the field can lead to disk instabilities which can increase 
the effective value of a within the disk. As an example, we can consider the MSI 
due to an external magnetic field. As we have already seen, such a field can give rise 



to a turbulent viscosity ~ VaH, where Va is the Alfven velocity due to the external 
field only. When (Va/c s ) > (H/r) 4 ^ 3 this will be the dominant source of angular 
momentum transport. The MSI is an incompressible instability, which operates only 
when the magnetic pressure is smaller than the fluid pressure, but similar effects can 
arise from compressional modes (Kaisig, Tajima, and Lovelace 1992; Tagger, Pellat, 
and Coroniti 1992; Tagger, Henriksen, Pellat, and Sygnet 1990) for an arbitrarily 
cold disk, i.e. a low (3. In this case the effects of the disk fluid pressure become 
small and on dimensional grounds we expect the the effective turbulent viscosity 
to scale as V^/fi. We see that the crossover point between the two mechanisms is, 
as expected, roughly where the magnetic field pressure from the external field and 
the fluid pressure within the disk are comparable. Both mechanisms are strongly 
dependent on the strength of the external field. They may be important close to 
the central object, or close to the inner boundary of the disk, but they will quickly 
become unimportant as we move to larger radii. In the region where they are 
important these effects must compete with the tendency of the central object to 
exert a torque on the disk through the stretching of its magnetosphere (Ghosh and 
Lamb 1979; Mauche, Miller, Raymond and Lamb 1990). The torque exerted in this 
way is B z B e r/Y;. Assuming that B g ~ B z Q*/£l this will be a large, even overwhelming, 
effect whenever the vertical MSI or compressional mode effects are strong enough to 
compete with the turbulent transport generated within the disk by purely internal 
processes. 

5. Conclusions 

More than 20 years after the seminal theoretical paper by Shakura and Sun- 
yaev the transport of angular momentum in accretion disks remains a contentious 
topic. In spite of this several points have emerged that are likely to be important 
in any future theory. First, the direction of angular momentum transport is not 
guaranteed by appeals to turbulent mixing. This makes models based on turbu- 
lent convection distinctly unpromising. Second, magnetic shearing instabilities do 
guarantee an outward transport of angular momentum and are unavoidable given 
a magnetic field. We believe that this instability will play a critical role in future 
theories of ionized disks. The implications for cold disks are less clear. Third, this 
makes an understanding of disk dynamos absolutely critical. Fourth, the tendency 
for successful disk models to include scaling laws of the form a « (H/r) n argues that 
a purely local process, e.g. a dynamo driven by the MSI process or by magnetic 
buoyancy is unlikely to work. We expect that there must be a nonlocal trigger 
for the dynamo, e.g. global waves of some kind. Finally, if global waves are im- 
portant in any way then the interaction between a three dimensional disk and its 
environment has to be understood quantitatively. 

Here we have presented a model for a dynamo in accretion disks based on 
the action of internal waves propagating from the outer edge of the disk. As the 
magnetic field grows it gives rise to shearing instabilities which transport angular 
momentum outward and magnetic flux vertically. These instabilities will eventually 



saturate the dynamo when (3 is still quite large. The resulting angular momentum 
transport is consistent with observations, although the prediction is rather approx- 
imate and is not strictly equivalent to a local viscosity in any case. From the point 
of view of dynamo theory the interesting feature here is a dynamo model which is 
nonturbulent, in the sense that the turbulent motions present are not responsible 
for driving the large scale magnetic field. From the point of view of accretion disk 
theory the interest of this model is that it provides a way of understanding an- 
gular momentum transport in accretion disks without resort to phenomenological 
considerations. 

The authors are happy to acknowledge useful contributions from Min Huang, 
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supported, in part, by NASA through contract NAGW 2048. 
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